Polarization dependence of x-ray absorption spectra in Na^CoC^ 
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In order to shed light on the electronic structure of Na^CoCh, and motivated by recent Co L- 
edge X-ray absorption spectra (XAS) experiments with polarized light, we calculate the electronic 
spectrum of a C0O6 cluster including all interactions between 3d orbitals. We obtain the ground 
state for two electronic occupations in the cluster that correspond nominally to all O in the -2 
oxidation state, and Co +3 or Co +4 . Then, all excited states obtained by promotion of a Co 2p 
electron to a 3d electron, and the corresponding matrix elements are calculated. A fit of the 
observed experimental spectra is good and points out a large Co-O covalency and cubic crystal field 
effects, that result in low spin Co 3d configurations. Our results indicate that the effective hopping 
between different Co atoms plays a major role in determining the symmetry of the ground state in 
the lattice. Remaining quantitative discrepancies with the XAS experiments are expected to come 
from composition effects of itineracy in the ground and excited states. 
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INTRODUCTION 



In recent years, there has been intense research on 
doped cobaltates because of their interesting transport, 
magnetic and superconducting properties. Na^OoC^ 
shows an exceptionally high thermopower over the dop- 
ing range 0.5 < x < 0.9, while also displaying low 
resistivity and low thermal conductivity 0, @- These 
properties are those wished in materials with poten- 
tial technological applications in refrigeration. In 2003 
Takada and coworkers found a superconducting state in 
Nao.3Co02 ■ 1.3H 2 at temperatures lower than 5 K 
The system Na^CoC^ is also interesting because of its 
similarities to the high-T c copper oxide superconductors, 
since both systems contain alternating layers consisting 
of oxygen and spin- 1/2 transition metal ions separated 
by layers of lower conductivity in an anisotropic crys- 
tal structure. Na^OoC^ consists of alternating layers of 
C0O2 and sodium, and becomes superconducting only 
when hydrated with water (Nao.35Co02 ■ l.SHgO)^]. 
However, in contrast to most cuprates, the C0O2 layers 
consist of edge-sharing CoOg octahedra, arranged in such 
a way that the cobalt and oxygen ions form a triangu- 
lar lattice. The superconducting resonance-valence-bond 
state for which some indications exist in the square lat- 
tice of the cuprates 0, , is believed to be favored in a 
triangular lattice 

In the most simple picture, the effect of sodium dop- 
ing is only to modify the nominal ratio of Co 3+ to Co 4+ 
depending on the sodium content. This in turn, as in 
the cuprates, modifies the number of free carriers and 
affects dramatically the conducting and magnetic prop- 
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erties of the system, x — corresponds to the highly cor- 
related (charge transfer or Mott) insulating limit. Mag- 
netic susceptibility measurements 0, have shown evi- 
dence of antifcrromagnetic long-range order below 20 K 
in a doping range of 0.75 < x < 0.9. From anisotropic 
dc magnetic susceptibility and /iSR measurements on a 
Nao.82Co02 crystal one concludes that the Co spin is ori- 
ented along the crystal c axis 0. 

Detailed knowledge of the electronic structure of this 
system is of course important for the understanding of 
all its properties. Ideally, from the knowledge of im- 
portant parameters such as the d — d Coulomb inter- 
action which includes exchange, the different hopping 
terms, and the charge-transfer energy one directly gains 
information about the band gap and the character of the 
electronic structure as it is explained in the classification 
scheme done by Zaanen, Sawatzky, and Allen |8j- Addi- 
tionally, knowing important orbitals and interactions in 
an intermediate energy scale accessible in optical exper- 
iments, like the three-band Hubbard model in the case 
of cuprates Q , one would derive an effective low-energy 
model (like the t — J model based on Zhang- Rice singlets 
[iol ITU IT2I IT3L ITU for the cuprates) which describes the 
relevant low-energy physics. As an example, interesting 
phenomena have been explained by the formation of the 
Zhang-Rice singlets, such as the unusual magnetic and 
transport properties in (La, Ca) x Sri4_ x CU24O41 M 

and 

other cuprates. Similar low-energy reduction procedures 
have been followed for the nickelates [16( and magnetic 
double perovskites ^3 • Instead, a systematic low-energy 
reduction has not been carried out for Na^OoC^ and it 
seems that no consensus exists on the low-energy effec- 
tive theory, in spite of the experimental and theoretical 
effort concerning the electronic structure of the system. 
At energy scales of 1 eV or larger, valuable information 
on this electronic structure is provided by local spectro- 
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scopic probes, like x-ray absorption spectroscopy (XAS) 
|l8tll9j . and x-ray photoemission (XPS) po|. while angle 
resolved photoemission (ARPES) [2ll l22l l2,'"il| brings in- 
formation on the dispersion relation just below the Fermi 
energy. In particular XAS, including a recent study of 
its dependence on the polarization of the incident light 
[l8L fl9| . indicate that Na^CoC^ exhibits the character 
of a doped charge-transfer insulator [8j, rather than a 
doped Mott-Hubbard insulator, and the ground state 
contains a significant O 2p portion. Several calculations 
of the band structure use a single band effective model 
H II |H HI H3 . These descriptions can be justified if 
the lattice distortion or crystal field break the degeneracy 
of the 3d ti g orbitals which are otherwise degenerate in a 
cubic environment, and if the O degrees of freedom can be 
integrated out in a similar way as in the cuprates. How- 
ever, Koshibae and Maekawa based on an estimation of 
the effects of the above mentioned lattice distortion, and 
simple geometrical arguments, proposed an alternative 
scenario, with four interpenetrating Kagome sublattices 
hidden in the C0O2 layer j2^. This scenario was further 
supported by Indergand et al. who derived the different 
hopping elements and interactions in the ensuing multi- 
band model |29j. Starting from this picture, Khalliulin et 
al. argued that the spin-orbit coupling of the correlated 
electrons on the ti g level, although relatively small (~ 80 
meV jU) strongly influences the coherent part of the 
wave functions and the Fermi surface topology at low 
doping and favors triplet superconductivity, in marked 
contrast to the cuprates |3fij . 

In this work, we solve exactly a cluster composed of a 
Co atom and its six nearest-neighbor O atoms for the ge- 
ometry and oxidations states relevant to Na^CoO^. The 
aim is to learn more about the electronic structure of 
the system in an intermediate energy scale by compari- 
son with experimental results on the polarization depen- 
dence of XAS spectra 0] . We confirm previous results 
concerning the highly covalent Co-0 bond and the low 
spin state of Co [ljj . In addition an analysis of the 
dependence of the XAS intensities on the polarization 
direction of the incident light, support the scenario of 
Koshibae and Maekawa that the effective Co-Co hopping 
through the intermediate O sites is an essential part of 
the low energy physics. Additionally it will be shown, 
that due to the strong covalency a large amount of holes 
reside at the oxygen sites. 

The formalism and relevant equations are presented in 
Section II. The results and its comparison to experiment 
are included in Section III. Section IV contains a sum- 
mary and discussion. 



II. METHOD 

A. The model 

An adequate starting point to describe the electronic 
structure of Na2,Co02 is a multiband Hubbard model in- 



cluding all Co 3d orbitals and O 2p orbitals, and their 
interactions (like the three-band Hubbard model for the 
cuprates @)- This model cannot be solved exactly and 
one has to resort to some approximation. Calculations 
starting from first principles approximate the interactions 
and are known to lead to wrong results when the on- 
site Coulomb repulsion Ua between transition- metal 3d 
electrons is important. For example, they predict that 
LaCu04 is a paramagnetic metal, while experimentally 
it is an antiferromagnetic insulator. In this work we take 
the opposite approach and retain exactly all 3d-3d inter- 
actions at the price to consider only one Co atom and its 
six nearest- neighbor O atoms, with octahedral Oh sym- 
metry. The effects of the distortions are discussed in sec- 
tion III. We also make some additional approximations 
which are not essential but reduce the size of the matri- 
ces of the Hamiltonian in each symmetry sector. This 
allows us to have a calculational tool flexible enough to 
allow a fit of the XAS spectra in a reasonable amount of 
computer time. For example, wc neglect the interactions 
between O 2p holes and this allows us to take us advan- 
tage of a restricted basis for these orbitals which has the 
same symmetry as the 3d tig and e g orbitals. The other 
O 2p orbitals become irrelevant in the absence of 0-0 
interactions. For the description of the Co L-edge XAS 
spectra, we also need to include in the Hamiltonian the 
energy of the Co 2p core hole and its repulsion with the 
Co 3d holes H mil. 

The Hamiltonian that describes the dynamics of holes 
inside the Co0 6 cluster is 

H = y ] 6t 3g d^ a d acr + ^ ' 6 eg dg a d at7 + ^""^ e j c jm C jm 

+ COPfaaPiav + ^ ^ (Pia<7 d 0<r + H.C.) 
iota i^jaf3a 

+u cd 

jmaa 

Here pf aa creates a hole on the O 2p orbital a at site i 
with spin a. The operator d^ a has an analogous meaning 
for the Co 3d orbitals. Similarly Cj m creates a Co 2p core 
hole with angular momentum j and projection m. 

The first two terms of Eq. describe the energy 
of the tig and e g Co 3d orbitals split by a crystal field 
e *2 S — £ e g — 10.Dg. We note that this difference can be 
called the "ionic" contribution to the crystal field split- 
ting, since after hybridization with the O orbitals, the dif- 
ference between mixed ti g and e g orbitals is much larger 
due to what is called the ligand field contributions. The 
third term describes the energy of the Co core hole, with a 
splitting £3/2 — 61/2 ~ 15 eV. The fourth term corresponds 
to the energy of the O 2p orbitals. The next two terms 
represent the Co-0 and 0-0 hopping, parametrized as 
usual, in terms of the Slater-Koster parameters 0| ■ The 
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term before the last includes all interactions between 3d 
orbitals originated by the Coulomb repulsion of electrons 
in the 3d shell. The last term describes the repulsion be- 
tween the Co 3d holes and a Co 2p core hole if present. 
We have neglected the exchange and higher multipole 
Coulomb interactions between these orbitals. Their ef- 
fects will be discussed later. Spin-orbit coupling of the 3d 
electrons is negligible compared with the other energies 
in the problem (~ 80 meV [3l]) and was also neglected. 

The matrix elements of the interactions inside the 3d 
shell are given by |35| 



Vx fiV p = / dr 1 dr 2 ipx(ri)(pfj,(r 2 )-, :(p„(ri)ip p (r 2 ), 

J |ri — r 2 | 

(2) 

where ip\(ji) is the wave function of the spin-orbital A. 
We have calculated all these integrals in terms of the 
Slater parameters Fn , F b , and F4 using known methods 
of atomic physics ti3 • The resulting form of this in- 
teraction term is long and we do not reproduce it here. 
The terms including only t 2g orbitals are described in Ref 
[38j |. and those among e g ones are described for example 
in Re f. I l6l | . Some of the remaining terms are listed in 
Ref. 37]. In particular, the value of the Coulomb re- 
pulsion between electrons of opposite spin at the same 
orbital is U c i — F + 4F 2 + 36F 4 , and the spin-spin inter- 
action between 3d orbitals of the same irreducible repre- 
sentation is J t . 2g = 3F 2 + 2OF4, and J Cg = 4F 2 + 15F 4 . 

In Na x Co0 2 , the formal oxidation state of the O ions 
is -2, and that of Co depends on the sodium content x, 
being Co 3+ for x = 1, and Co 4+ for x — 0. This means 
that our cluster has to be solved for a total of four holes 
(configurations d 6 L°, d 1 L 1 , d 8 L 2 , etc.) in the former 
case and five holes (configurations d 5 L°, dPL 1 , etc.) in 
the latter. The respective probabilities are assumed to be 
x and 1 — x. For simplicity, we will label these two cases 
as Co 3+ and Co 4+ respectively, although as we shall see 
there is a significant degree of covalency and the real Co 
oxidation states are smaller. 



leads to the operator 

H L = J2(0+ z<7 d 2lza + R.c.). (3) 

Here the destruction operators of the 3d holes d 2 i z<7 are 
expressed in terms of Z z , the orbital angular momentum 
projection. is an irreducible operator that trans- 

forms like an angular momentum L = 2 with projection 
l z , constructed combining the components of A in spher- 
ical harmonics (Au z ) with the core hole operators a 
that create a Co 2p hole with angular momentum projec- 
tion l z and spin a . The expression of the five components 
of O21 a {l z = —2 to 2) is derived easily using Clebsch- 
Gordan coefficients. Similarly, using these coefficients, 
it is easy although lengthy, to express the light opera- 
tor Hl in terms of the basis for the orbitals used in the 
Hamiltonian (Eq. (Q) and Cartesian coordinates of A: 

'Hl = ^ a PpnaaAi3P^ m d aa + H.c. (4) 
ftjmaa 

This is the most convenient form for our purposes. 

Using Fermi's golden rule and neglecting an unimpor- 
tant prefactor, the XAS intensity becomes 

/ = ^p i |(/|W i |i)| 2 ^ / -^-^). (5) 
if 

Here \i), is one of the two possible initial states (ground 
state for Co 3+ or Co 4+ ), pi its probability that depends 
on the Na content, and Ei its energy. Similarly / labels 
the final states, fku is the energy of the incoming light 
represented by the vector potential A. 

In order to simulate the measured spectra, we have 
broadened the delta functions appearing in Eq. © 
replacing them by a Lorentzian line shape with a full 
width of 0.5 eV at half maximum (FWHM) at the 
£3-edge and 0.7 eV at the i 2 -edge due to different life 
time broadening effects. 



B. The x-ray absorption intensity 

The ground state of the cluster in both above men- 
tioned cases (for Co 3+ and Co 4+ ) does not contain any 
Co 2p core hole. The excitation with light promotes a 
Co 2p electron to the Co 3d shell, or in other words, 
creates a core hole and destroys a 3d hole. In the dipo- 
lar approximation, the effect of the light is calculated in 
time-dependent perturbation theory from the addition of 
the term Hl = P ■ A to the Hamiltonian, where p is the 
momentum operator and A = (A x , A y ,A Z ) is the vec- 
tor potential. Except for an unimportant prefactor, the 
relevant part of TLl can be deduced from symmetry ar- 
guments and general physical considerations: it should 
be rotationally invariant, spin independent, and contain 
terms that destroy a 3d hole and create a core hole. This 



III. RESULTS 

A. The structure of the ground state 

There are already indications that the Co ions in 
Naa;Co0 2 are in a low spin state |40L l4ll| . supported also 
by the interpretation of XAS measurements using polar- 
ized light [l8L TTg| . This means that crystal field effects 
dominate over the exchange terms. In fact, from the pa- 
rameters that best fit the experiment (see Table lHl|) . one 
sees that \QDq = 1.2 eV. As mentioned in the previous 
section, this corresponds to the "ionic" part of the crystal 
field. The actual splitting is much larger due to the cova- 
lency effects which are also strong in the system. In table 
[1 we present an overview of the configurations present in 
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Co 3 + 

Groundstate 


Co 4 + 
Groundstate 


d 6 


0.30 


d 5 


0.14 


d 7 L* 


0.47 


d 6 L J 


0.40 


d 8 L 2 


0.20 


d 7 L 2 


0.35 


d 9 L 3 


0.03 


d 8 L 3 


0.10 


d 10 L 4 


0.00 


d 9 L 4 


0.01 






d 10 L 5 


0.00 



TABLE I: Overview of the hole distribution for the ground 
states. The numbers give the hole distribution in percent. 



the ground state for clusters with nominal oxidation state 
Co 3+ and Co 4+ . 

The ground state for Co 3+ is a singlet with A\ g symme- 
try (for the sake of clarity we use capital letters to denote 
the irreducible representations of many-body states). For 
the cluster representing nominal Co 3+ , the most impor- 
tant configurations in the ground state are d 6 (30%) and 
d 1 L 1 (47%), indicating a very strong covalency in the 
system. Note that the d 6 population is smaller than the 
d 7 L population and considerably less than 50% of the 
total. Such a behavior can only be explained if one takes 
more than one ligand hole into account. In a configura- 
tion containing only d G and d 7 L with a positive charge 
transfer gap A CT = E{d n+1 L) - E{d n ) (see table III), 
the d e population will always be larger or equal to d 7 L, 
but will never be smaller. Things change if one includes 
also configurations containing two ligand holes d 8 L 2 with 
sufficiently large hopping terms. In this case the relative 
hole population changes towards the d 7 L configuration, 
so that even for a positive charge transfer gap Act the 
dominant configuration can include a ligand hole. This 
is a rather important result since it demonstrates that 
the hybridization between all those states involving one 
or more ligand holes is large enough to produce a low 
energy state which in fact is lower than that of the start- 
ing d 6 configuration containing no ligand holes. The d 6 
configuration corresponds to t\ g (all holes with e g sym- 
metry occupied). For simplicity to reduce the size of 
the matrices in the fitting procedure, we neglect the con- 
figurations ^2g e g- This is justified by the fact that the 

6 

■2<J 



interaction term which mixes the 3d configurations $ 

and t\ g e\ is A = V3(F 2 - 5F 4 ) ~ 0.3 eV (see Table ITTH. 
is much smaller than the effective crystal-field splitting 
~ 3 eV. The amount of t\ g e\ 3d configurations can be 
calculated by perturbation theory, but in any case, its 
influence on the XAS spectra is of the order of 1%. Sim- 
ilarly, for Co 4+ , we have neglected states with more than 
one 3d t^g hole. 

With these simplifications and because of the neglect 
of interactions between electrons at the O sites, the num- 
ber of relevant states in the subspace in which the ground 
state for Co 3+ is, reduced to 7 singlets with A\ g sym- 
metry. Similarly, the ground state for Co 4+ is obtained 
from six identical 38x38 matrices and corresponds to a 



d6 t 2q5g q l L e q l 



(a) 



(b) 



C ' 7 t2g 6 eg 1 '-eg'-t2 





d VLt 



:• ; 



d5 t2g 5 e g l'-eg 



(c) 



FIG. 1: Sketches of the energy scheme for different hole con- 
figurations of Co 4+ . (a): Hybridization between states with 
none and one ligand hole, (b) and (c) : Hybridization between 
states with one and two ligand holes. 



spin doublet with Ti g symmetry (this means three-fold 
orbital degeneracy). 

Due to the symmetry of the cluster (Oh), the 3d holes 
conserve its symmetry when they hop to the O 2p states. 
Neglecting the small amount of 3d i| s e g states in Co 3+ , 
the only possibility for a ligand hole is then a linear 
combination with e g symmetry giving the configuration 

t%„el el 



Co 4+ (five holes) 
Ground state 



d 5 


0.14 


d 6 L c \ 


0.37 


d?L k 


0.28 






t 2g 


0.04 


V2g 


0.07 



TABLE II: Hole distribution for the ground state of the Co 4+ 
cluster, including information of symmetry up to configura- 
tions with two ligand holes. The numbers give the hole dis- 
tribution in percent. 

A more complicated picture is expected for the cluster 
corresponding to Co 4+ , since it contains five holes mean- 
ing that there will always be one hole with ti g symmetry 
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in the system (either at the Co site or distributed among 
the O ions). For the parameters that fit the experiment, 
the ground state of this cluster is a spin doublet with 
orbital symmetry Ti g . This three-fold degeneracy in the 
point group Oh is split by the rhombohedral distortion 
as discussed in the next subsection. While an overview 
of the structure of the ground state is given in table [I] 
the distribution among different possible configurations 
up to two ligand holes for the ground state of the Co 4+ 
octahedra is shown in table |H] Note that states with 
one and two ligand holes dominate the composition of 
the ground state. The d 5 configuration is very simple. It 
consists of four 3d eg holes and one 3dt 2g hole, or equiva- 
lently it consists of five 3d t2a electrons (d% ). The con- 

9 Z 2g 

figurations with one ligand hole with the same symme- 
try are d% 1 L e i, and d% L f i . The former has a larger 

l 2g e g 9 l 2g 2 9 

hybridization with the d 5 configuration because the hop- 
ping between 3d and 2p e g or ti g electrons involves the 
Slater-Koster parameter (pda) or (pdn) respectively, and 
the former has a larger magnitude. For configurations 
with two ligand holes there are two configurations with 
Ti„ symmetry: dj 6 x L e \L t \ and dj 5 2 L e 2. They both 

t 2g e g 9 2 S t 2g e g S 

hybridize with the dominant one hole state d% 1 L e i. 

2g e g 9 

The only difference is that the first configuration couples 
via tt 2 hopping to the one hole state while the second 
configuration couples via a t e hopping which is larger 
and results in a stronger occupation of this state (see ta- 
ble [HJ|. Although it is allowed by symmetry, the above 
mentioned states with two ligand holes do not hybridize 
with the state d% L t „ due to the particular structure of 

l 2g 

the wave functions. 

^From the above described distribution of holes in the 
Co 3 and Co 4 clusters, it turns out that the real Co valence 
for each cluster is 2.04 and 2.56 respectively, emphasiz- 
ing again the very strong degree of covalency in these 
systems. 



B. Breaking of the octahedral symmetry 

The edge-sharing CoOg octahedra in a C0O2 layer of 
Na x Co02 are compressed along the c axis (oriented par- 
allel to the (111) direction of one octahedra, see Fig. 
2(a) I. This distortion reduces the point group sym- 
metry of our cluster to D^d- As explained above, the 



ground state of the Co 4+ cluster belongs to the T 2g rep- 
resentation of Oh- This representation is split in D$d as 
Ti g = A' lg + E' g , where to avoid confusion, the represen- 
tations of D%d are primed. Labelling as \xy), \yz), and 
\zx) the basis functions of the representation T 2g (which 
in our case correspond to the many-body wave functions 
of the Co 4+ cluster), the basis functions that transform 
irreducibly in D^d can be written as 




FIG. 2: (a): Distortion of the octahedra. (b): Splitting of the 
ti g states in a trigonal crystal field. 



for the state invariant under operations in D$d and 

\E' g ±) = -L(\xy) + \yz) + e ±i * \zx)) (7) 

for the doubly degenerate states that transform like E' . 
A 3d orbital with the symmetry A' lg looks like a 3z 2 — r 2 
orbital and points into the direction of the crystallo- 
graphic c axis. In a one-electron picture, for an adequate 
ordering of the energy of the orbitals, the Co 4+ cluster 
would correspond to e' g orbitals completely filled with 4 
electrons and a half filled a' lg orbital ■ 

The rhombohedral distortion of the C0O6 octahedra is 
estimated by the deviation of the Co-O-Co bond angle 
from 90°. For example for x w 0.7 Huang et al. find 
an angle a ~ 96° |46j which leads to a change of the 
thickness of one octahedra by 10%, assuming that the Co- 
O bond distance remains constant during the distortion. 

/From the fit of unpolarizcd XAS spectra [18j . Wu 



et al. estimated D 



try 



1 eV for the splitting D 



try 



1 



y/3 



(\xy) + \yz) + \zx)) 



(6) 



(E' g ± \H\E' g ±) - (A' lg \H\A' lg ) between the many-body 
ground state A' lg and the excited states \E' g ±). Instead, 
Koshibae and Maekawa using a point charge model ob- 
tained a negligigle splitting and suggested that the band 
splitting at the zone center obtained in band structure 
calculations is dominated by the effective Co-Co hopping 
mediated by intermediate O atoms 28]. However, one 
expects that the splitting is dominated by covalency ef- 
fects and that the ionic contribution plays a minor role. 
Very recently D trg = 0.315eV has been obtained by 
an ab-initio quantum chemical configuration interaction 
method in a CoO@ cluster |44|. This is the most reliable 
result. The method used is known to provide accurate 
results for charge conserving excitations and has been 
used succesfully in a large family of strongly correlated 
systems [44^. 

To check if the distortion affects our results for cu- 
bic symmetry, we have estimated the matrix element 
Q mix = (Eg ± \H\E' g ±,2) between the E' g states and 
the lowest exited states \E' g ±,2) of the same symme- 
try. These excited states belong to the irreducible repre- 
sentation E g of Oh and therefore do not hybridize with 
the ground state in cubic symmetry. The only part of 
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the Hamiltonian where the trigonal distortion comes into 
play is a change in the hopping part through changes in 
relative orientations of the orbitals, and 0-0 distances. 
We assumed an r~ 3 distance dependence for the hop- 
ping between p orbitals [4j|. We find Q m ix — —0.09 eV. 
Since |Q m ia;| is very small compared to the other ener- 
gies of the model, in particular IQDq and the separation 
between Ti g and E g levels, we can safely use the states 
obtained in cubic symmetry for the calculation of the 
XAS spectra. 

A ground state with E' g symmetry is incompatible with 
the XAS data. Therefore only a positive D trg is con- 
sistent with experiment. The value D trg — 0.315eV is 
smaller but of the order of magnitude of the band width 
reported in ab-initio calculations. A small value of Dt rg 
would be consistent with the proposal of Koshibae and 
Maekawa that the energy splitting does not originate 
in the crystal field, but is determined by the kinetic en- 
ergy of the electrons. Within this picture, translated to 
our many-body states, and neglecting D trgi the ground 
state for one Co 4+ cluster moving in a lattice of Co 3+ is a 
Bloch state with total wave vector K — in an effective 
Kagome sublattice with point symmetry A lg . To keep it 
simple, restricting the wave function to only three sites, 
and one of the four Kagome sublattices, the ground state 
would be 



\g) = -j=(\l,xy) + \2,yz) + \3,zx)), (8) 

where here \i,a) denotes a state composed of the ground 
state of one Co 4+ cluster with symmetry a at site i, and 
the (invariant A\ g ) ground state of the Co 3+ cluster at 
the other two sites. However, the effective tight-binding 
Hamiltonian considering only the most important hop- 
ping term and neglecting D trg , has a U(4) symmetry due 
to the equivalence between the four Kagome sublattices, 
and the different possibilities of breaking of this symme- 
try leads to several scenarios for the ground state [2j| . 

The difference between Eqs. © and JSJ) for the polar- 
ization dependence of the intensity in the XAS spectra, 
can be easily illustrated for the case in which the final 
states \f a ) in Eq. (JSJ for a Co 4+ cluster have a core hole 
2p a (a — x, y, or z) with given spin, and the rest of 
the state (the holes of the 3d shell) belongs to the A\ g 
irreducible representation. This intensity with a factor 
two coming from spin, corresponds to the sum of the in- 
tensities of the peak A in Fig. 3 (for a total momentum 
j = 3/2 of the core hole) and the corresponding peak 
at higher energies and lower intensity for j = 1/2 (A' 
in Fig. 3). For incoming light in the direction (3, TLl is 
proportional to the momentum operator in this direction 
pp. Then by symmetry, all matrix elements en- 
tering Eq. (0 can be related to 7 = (fx\p y \xy). Using 
Eq. ijfjj , it is easy to see that the contribution of these 
states (summed over a) to the intensity in the one-cluster 
picture, for A is parallel and perpendicular to the c axis 



becomes 

4=4| 7 | 2 /3,/I = | 7 | 2 /3. (9) 

In the extended picture, naturally, the light operator is a 
sum over all sites of that already described. Taking again 
only three sites for simplicity, there are three possible 
final states for each core hole orbital 2p Q \i, f a ) depending 
on which site i was excited by the core hole. Adding these 
contribution one finds, using Eq. (JSJ) 

jext = jext = 2| 7 | 2 /3. (10) 

Thus due to the loss of coherence, there is a factor 2 
less intensity for A parallel to c in the XAS spectra for 
the sum of the contributions of peaks A and A' (see Fig. 
3) in the case in which the hole is distributed over one 
Kagome sublattice. This factor is large enough to be de- 
tected by the experiment. In the latter case, the result 
is the same as for a ground state composed of statistical 
average of the three T^ g states and has no polarization 
dependence. The situation is different if the breaking of 
the above mentioned U(4) symmetry leads to some co- 
herence between Ti g states at the same site. Note that 
the intensity is redistributed among the different direc- 
tions but keeps the same angular average. This is due 
to the local character of the transition and the fact that 
only one irreducible representation involved in the initial 
state, final state and operator. 

We shall show that our results point to an intermediate 
situation, but with dominance of the lattice effects over 
the on-site distortion. 



C. The XAS intensity 

As mentioned above, the ground state for a Co 3+ clus- 
ter is invariant under the symmetry operations of Oh 
and as a consequence, the corresponding XAS intensity 
is independent of the polarization of the incident light. 
Also, for Co 4+ , the XAS intensity for a statistical aver- 
age of the three states of the irreducible representation 
T2 g {\xy), \yz) and \zx)) (which would correspond to a 
ground state like Eq. (JSJ) leads to the same intensity for 
A||c and A_L c. 

Instead, using Eq. J§J| for the ground state of Co 4+ , 
one obtains a polarization dependence similar to the 
experiment, but more pronounced: The ratio of the 
intensities for excitations to a state with A\ g symmetry 
plus a core hole is — 4 (sec Eq. ©) whereas the 

experimental ratio has been found to be I^ xp /I e * p w 1.8. 
Therefore the experiment suggests an intermediate 
situation between the ground states © and (JSJ. The 
actual many-body ground state of the system is out 
of the scope of our cluster calculations. In order to 
simulate this situation, we have calculated the spectrum 
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FIG. 3: XAS spectra for different polarizations: (a) E\\c and 
(b) E _L c Dashed line: Co 3+ cluster. Dotted line: Co 4+ 
cluster. Straight line: sum of both representing (Nao.sCoO^). 
In (a) and (b) Eq. © was used for the ground state of the 
Co 4+ cluster was used, while (c) corresponds to Eq. JSJ giving 
a polarization independent result. 

as a weighted average of Eqs. © and i|lUfl: 



Itotal = (l-Olf+ZIf Xt (11) 

with < £ < 1. In order to get the observed experimental 
ratio for the Ai g peak £ has to be approximately £ = 0.73. 

In Fig. |3 we show the results for the XAS intensity 
for the parameters indicated in Table II I II that best fit 
the experiment. As adjustable parameters, we have used 
F (or U d = F + AF 2 + 36F 4 ), e Q (or A CT (Co n+ ) = 



F = 3.5 pda = 2.35 WDq i = 1.2 

F 2 = 0.2 pd-K = -1.0 e p = 13 

F 4 = 0.006 ppa = 0.8 A C t(Co 3+ ) = 2.9 

U dd = 4.52 ppn = -0.2 A C t(Co 4+ ) = -1.1 



TABLE III: Best fit parameters, all parameters given in eV. 

E(d n+1 L) - E{d n )), lODq ,(pda), and (ppa). We have 
taken Udc = Ud + 0.15 eV. F 2 and F4 were taken from 
values that fit approximately the energy separation be- 
tween different terms in atomic spectra |3!j. We as- 
sumed the relations (pdir) = — \^3(pda) / 4, and (ppir) — 
-0.81(pdcr)/3.24 |H. 

In terms of more fundamental parameters one has: 

A(Co 3+ ) = eo - e eg - 3^0 + 8^2 - 33F 4 - ppa + ppir 
A(Co 4 +) = e Q - e eg - 4F + AF 2 - 37F 4 /2 - ppa + ppir, 

(12) 

and correspond to the minimum energy necessary to pro- 
mote a Co 3d hole to a linear combination of O 2p states 
in absence of Co-0 hopping. If the O 2p hole is localized 
at one site, ppa and ppn are absent in these equations 
and both A(Co 3+ ) increase in 1.0 eV. 

To fit the experiments and explain the different ob- 
served structures, several configurations with consider- 
able amount of ligand holes (described above) should be 
present in the ground state. This implies large covalcncy, 
and this in turn points to high hybridization, a low charge 
transfer gap Act and low Ud in comparison with other 
transition metal oxides. Also, the relative position of the 
structures for Co 3+ and Co 4+ imply relatively large lODq 
and low spin configurations. This is a rough explanation 
of the resulting parameters in Table IIIII 

The Co 3+ cluster (dashed line in Fig. |3J) shows only 
one main peak at the L^^ edge (peak (B) in the figure 
EJ) which results from an excitation from the ground 
state with A± g symmetry to an excited state with E g 
symmetry and is polarization independent. Comparing 
to the similar system LiCo02 which nominally contains 
only octahedra with four holes (i.e. a Co 3+ -central ion) 
|47j|. one observes a weak structure on the high energy 
side of the e g peak which is missing here because of 
the missing 2p-3<i core hole exchange interaction. The 
Co 4+ cluster shows two main structures at both edges. 
The lower energy peak structure (A) contains of a single 
peak and results from excitations into states with Ai g 
symmetry, while the higher energy peak structure (C) 
originates from many excitations of low intensity to final 
states with Ti g and T 2g symmetry. Only the former 
peak (A) shows a polarization dependence and is strong 
for incoming light polarized along the (111) direction 
(crystallographic c axis) and weak for polarization per- 
pendicular to the (111) direction (in- plane polarization). 
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FIG. 4: Comparison of the theoretical results (bottom) with experimental data (top) for x=0.4 and x=0.6 for polarizations 
indicated inside each figure. In the middle curve the results are manipulated in such a way, that the intensity of the Co 4+ 
cluster is decreased by a factor 2.0. 



In Fig. 0]we compare the theoretical results with ex- 
perimental data for the stoichiometries Nao.4Co02 and 
Nao.6Co02 as well as for polarization perpendicular and 
parallel to the crystallographic c axis. In the theoretical 
curve at the bottom, we have used the weighted average 
Eq. with £ = 0.73 which best fits the polarization 

dependence of peaks A and A'. As explained above, this 
average redistributes the spectral weight among the dif- 
ferent polarizations but does not affect the total intensity 
coming from a Co 4+ cluster. As it is apparent comparing 
the bottom and top curves in Fig. ^ the above mentioned 
average is not enough to reach a quantitative agreement 
with experiment: there is still a disagreement in the in- 
tensity ratio between that coming from a Co 4+ cluster 
and that originating from a Co 3+ cluster. The theoret- 
ical contribution of Co 4+ seems to be two times larger 



than in the experiment for x=0.4 (this ratio is a little bit 
less than two for x=0.6). Wet chemical redox analysis 
point to a proportion of Co + smaller than that corre- 
sponding the nominal Na content |50| . This can explain 
part of the discrepancy. Another possible explanation for 
the discrepancy is the redistribution of spectral weight 
due to itineracy of the final states. The effects of the 
hopping in a model of two sites has been studied by GA. 
Sawatzky and A. Lenselink |4^. In our case, these effects 
seem to be of minor importance for the peak A, because 
the expected magnitude of the hopping ~ 0.1 eV [2^ is 
much smaller than the energy difference between excited 
states involved (~ 1.9 eV). However they are certainly 
important for the higher energy multiplets. 

Note that our neglect of the exchange interaction be- 
tween the core hole and the 3d holes leads to slight change 
in the shape of the curves, but does not affect the total 
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intensities for both clusters. 

Spin-orbit coupling inside the 3d shell, which we have 
neglected might also lead to a reduction of the intensity of 
the low energy peak at the L^-edge that corresponds to 
peak A as an effect of selection rules 0|- However, in our 
case, the effect of spin-orbit coupling is small compared 
to the intersite hopping terms t ~ 0.1 eV [2^. Therefore, 
our approach of considering a sixfold degenerate ground 
state for Co 4 (T eg doublet) instead of split doublet and 
quartet states is justified. Therefore, we expect that the 
reduction of intensity due to spin-orbit coupling inside 
the 3d shell is small. 

The relative peak positions do not change upon dop- 
ing, except for the Ai g peak. The distance between the 
main (E g ) peak and the A± g peak becomes smaller with 
increasing x (decreasing amount of Co 4+ ions) 19]. This 
change is best simulated by a variation of the ionic crys- 
tal field lODq from 1.2 eV (x=0A) to 0.9 eV 0=0.6). A 
decreasing value of lODq is also supported by Huang et 
al. who performed neutron diffraction on powder sam- 
ples and found an increasing Co-0 bond length with in- 
creasing doping . This would change the total crystal 
field splitting (the ionic term lODq plus the effects of 
hopping), but a slight change in the hopping parameters 
does not change the relative peak positions. 

IV. SUMMARY AND DISCUSSION 

We have investigated the local electronic structure 
around a Co atom in Na x Co02, by solving exactly a 
Co06 cluster (a basic octahedron of the system) contain- 
ing the Co 3d and O 2p valence electrons, and with the 
appropriate charge according to the Na doping. We have 
included all interactions between 3d electrons in the clus- 
ter. For the calculation of the polarization dependence of 
the XAS spectrum, we have included the energy of the Co 
2p core hole with its spin-orbit coupling Aso ~ 15 eV, 
and its repulsion with the Co 3d holes. The exchange in- 
teraction between Co 2p and 3d holes has been neglected 
for simplicity. Due to comparatively large magnitude of 
Aso, crystal-field effects and covalency, one expect the 
influence of this term to be small in comparison to lighter 
3d transition metals [42| . 

Within a purely local picture, assuming an ground 



state [Eq. (JBJ] and disregarding for the moment the total 
intensity of the Co 4+ contribution, we are able to repro- 
duce the essential features of the polarization dependent 
XAS experiments fairly well. Slight differences at the 
shoulders of the largest peaks originate from the neglect 
of the Co 2p-3d exchange interaction which would add 
an additional weak structure in this region, as shown by 
comparison with the case LiCo02 |47j . 

Quantitatively, a fully coherent A' lg ground state, as 
Eq. jnj leads to too large intensity ratio for peak A 
(and A') of the Co +4 contribution between light polar- 
ized parallel or perpendicular to the trigonal axis. This 
points out that a purely local picture is not fully consis- 
tent with experiment. In fact, the best fit corresponds 
to a w 70% loss of coherence produced by delocalization. 
This redistribution of intensity among the different polar- 
izations is consistent with an itinerant picture based on 
four Kagome sublattices hidden in the C0O2 layer |2SL l29| 
as a first approximation for the electronic structure. 

The above mentioned redistribution is not enough to 
explain the large total amount of the Co +4 contribution 
to the intensity. Part of this discrepancy might be due 
to proportion of Co +4 smaller than that suggested by 
the nominal Na content [50|. A more quantitative cal- 
culation of the XAS spectrum requires to take into ac- 
count the itineracy of holes in the ground state and the 
excited states after the creation of the core hole, which 
leads to spectral weight transfer between Co +4 and Co +3 
contributions |49| . This is in principle possible if one 
combines the knowledge of the local transition matrix el- 
ements with an effective Hamiltonian for the motion of 
the holes through the lattice. This approach has been 
followed to calculate different optical properties of the 
cuprates.0, 0, |^, |U H3, 03 However, in this system 
this task is more difficult and is beyond the scope of the 
present paper. 
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